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We discuss an instability in a large class of models where dark energy is coupled to matter. In 
these models the mass of the scalar field is much larger than the expansion rate of the universe. We 
find models in which this instability is absent, and show that these models generically predict an 
apparent equation of state for dark energy smaller than -1, i.e., super-acceleration. These models 
have no acausal behavior or ghosts. 
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MOTIVATION 

Observations of distant Type la supernovae 0, 0| and 
the cosmic microwave background [3J] together strongly 
prefer an accelerated expansion of the universe in the re- 
cent past. In the standard cosmological model this is 
accommodated by introducing "dark energy" , a compo- 
nent which has a significantly negative pressure causing 
the expansion of the universe to accelerate. 

In the standard cosmological model, dark energy is 
completely decoupled from the rest of the matter in the 
universe except for its gravitational effects. It is interest- 
ing to consider more general models in which the dark 
matter and dark energy have a coupling. Such models 
could have new nontrivial signatures in cosmology and 
structure formation. 

One simple class of such models is a model in which the 
vacuum energy density depends on the matter density. 
We shall consider a class of these models in which the 
dark energy responds to changes in the matter density on 
a time scale shorter than the expansion time scale. For 
example, one can consider models with scalar field dark 
energy coupled to matter (e.g., 0, 0, @, 0, 0, H in 
which the mass of the scalar field is much larger than the 
expansion rate (for example, the MaVaN scenario 

As we show below, these models generically suffer from 
an instability which we label AZK-instability. The AZK- 
instability was pointed out in the context of mass- varying 
neutrinos (MaVaN) A similar effect was identified 

in the context of unified dark energy models (l3j . This 
instability can also occur in models of dark ene rgy cou- 
pled to matter, such as the MaVaN scenario [llj, the 
Chameleon dark energy scenario [3| and the Cardassian 
expansion scenario [lfj . Not all models in the above sce- 
narios are necessarily unstable (for example, ltl 17. 18|). 
This will become clear when we discuss the instability. 

In this paper, we will construct a large class of models 
in which this instability is avoided. We find that these 
models generically predict an apparent equation of state 
(pressure over energy density) wde which is less than -1 
(such a phase is labeled super-acceleration That 
is, a model of interacting dark energy can be incorrectly 
interpreted as a theory with super-acceleration if the in- 



teractions are not taken into account. 

For example, the coupling of dark energy to matter 
could be such that the total matter density decreases 
more slowly than 1/a 3 (where a is the scale factor of 
the universe). When we interpret observations in such 
a universe with a canonical matter density term (that 
decreases with expansion as 1/a 3 ) and dark energy, we 
would infer an equation of state for dark energy more 
negative than it truly is [2CJ, l2l|. There is no physical 
reason why this inferred equation of state cannot be be- 
low -1. 

This is particularly interesting because current data 
seem to favor a dark energy density which is almost con- 
stant or even increasing with time 22, 23, 24, 25, 2(| 27 . 
[2^. [29l . [30l[3]| . and exciting results can be expected in the 
future 1321133. 34, 35| . SNIa observations currently favor 
a phase of super-acceleration. Future SNIa and CMB ob- 
servations have the potential to detect super-acceleration 
[l9j |. No other combination has been shown to robustly 
detect the signature of super-acceleration, although com- 
bining SNIa and baryon oscillation [30| or weak lensing 
data set seem promising. Note that a measurement of 
just the average equation of state [3^ is not sufficient for 
this purpose 37 1. This was made explicit recently [38| 
using a simple single scalar field model. 

Scalar field models with canonical kinetic terms al- 
ways produce u>de > — 1- Effective models with the 
opposite sign kinetic term [H, [3t| imply wde < — 1 
but are unstable [i(| unless more than one scalar field 
[ill I42I [43l . H , [45( 1 or quantum effects [46| are considered. 
Models with higher derivative terms or scalar-tensor the- 
ories can give rise to a n ap parent wde < — 1 but 
are constrained 4^, 3j| 5(|. Interpreting an alternative 
gravity theory in the context of 4-d GR can also lead to 
super-acceleration [Hi, 52, 53, 54, HBl 5^1. Some Cardas- 



sian models may have wde < — 1 57, 53, 55^ while still 
satisfying the dominant energy condition. Another pos- 
sible way to get super-acceleration with no instabilities 
is to appeal to photon-axion mixing (conversion of pho- 
tons to axions) in a universe dominated by a cosmological 
constant (or quintessence) [fjoj . 

In our models, the superacceleration arises due to in- 
teractions of dark energy and matter. Our models there- 
fore provide super-acceleration with none of the atten- 
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dant problems that plague most of the above models. 
Furthermore, the interactions are generic; we do not need 
to fine-tune couplings in order to avoid theoretical pit- 
falls or observational constraints. We therefore believe 
that considering interactions of dark energy is the best 
way to generate models of superacceleration. 

AZK-INSTABILITY 

In this section we will consider a general class of models 
in which the dark energy density is coupled to the non- 
relativistic matter density. For an example of how this 
could occur, suppose that non-relativistic matter parti- 
cles are coupled to a scalar field. Thus the local density 
of the matter particles can influence the vacuum expec- 
tation value (vev) of the scalar field. The change in the 
potential of the scalar then affects the dark energy, thus 
coupling matter and dark energy. 

In this class of models, the matter fields will be taken 
to have a matter density tim- They are coupled to a scalar 
field x (dark energy) through Yukawa like couplings. We 
take the potential to be 

E = Jd 3 x V(x,n M ), (1) 
= J d 3 x [V {x) + mn M + \g(x)n M } ■ (2) 

We will assume that m 2 , = Vq(xo) + W(Xo) n M, the 
mass-squared of the scalar field about its vev xo, is very 
large so that the x field always sits at the minimum of its 
effective potential. This is the central assumption of our 
paper. The mass will certainly have to be larger than the 
expansion rate of the universe to be consistent with this 
assumption. We will also assume that the mass is large 
enough to satisfy the constraints imposed by experiments 
that probe the strength of a fifth force. 

In the absence of the last term, this is the potential en- 
ergy of two decoupled fluids. The first term corresponds 
to a cosmological constant term (since we have assumed 
that the field x is always at the minimum). The second 
term is the energy density of a dark matter fluid with 
density tt-m and particle mass m. 

The last term couples these two fluids, and leads to in- 
teresting effects. In particular xo, the value of the scalar 
field at its minimum is now found by solving the equation 

V^Xo) + Xg'(xo)n M -0, (3) 

where Vq and g' are derivatives of Vq and g with respect 
to x- Thus xo is now a function of tim- 

We can make the dependence of xo on um explicit in 
the following way. Consider small deviations in tim- The 
vev of the scalar field shifts to account for this change in 
riM- Taking a further derivative, we find 

(K'(xo) + W'(Xo)nu)^- + Xg'ixo) = . (4) 



This explicitly shows how xo varies as Hm varies. 

In writing Eq. [2l we neglected the kinetic term in com- 
parison to the potential. This is necessary if the scalar 
field is to behave as dark energy and, as we now show, 
consistent with our assumption of a large mass for the 
scalar field. Note that x = "-m^xo/^tim- Working out 
this expression, we find that x 2 /V for X — Xo is given 
by (VQ 2 /Vm^,)(n'M/m x nM) 2 - Lets look at changes to 
the scalar field potential around X = Xo- Unless there 
are strong fine-tunings and cancellations, we will have 
VqSx < V and m^(<5\;) 2 /2 < V, which together imply 
that 2Vq 2 /Vm 2 , < 1. Hence the natural expectation is 
that x 2 /y ~ H 2 /m 2 x . For large enough m x , the kinetic 
term is negligible. 

We now show that there is an instability in this system. 
We start with a configuration where the dark matter is 
evenly distributed, and the x field is at its minimum xo 
everywhere. Now consider small fluctuations in the mat- 
ter density Sum which preserve J d 3 x Sum — 0, i.e., the 
total number in volume r. The integral is over some re- 
gion r, much smaller than the Hubble volume, over which 
the fluctuations are coherent. Such a fluctuation leads to 
a change in the total energy. The energy change propor- 
tional to Sum vanishes because of Eq. [3] and the condition 
that J d 3 x Sum = 0. The energy change to next order is 

= -- / fc M A — (6) 

2 J m z x 

Therefore the leading correction to the energy is always 
negative, implying that the configuration is unstable to 
the growth of these fluctuations. We dub this the AZK- 
instability. This instability was first noted in the context 
of the MaVaN scenario [12j . 

We have neglected gravity and the expansion of the 
universe in the above analysis. We neglected gravity be- 
cause the relevant length scales are much smaller than 
the Jeans length; the instability occurs on all scales and 
hence the effect is most severe on microscopic scales. The 
analysis above was thus for a region r much smaller than 
that where gravity would be important. We neglected 
the expansion of the universe because the relevant time 
scales are much smaller than the age of the universe. In 
addition our setup started with a smooth distribution 
of matter. For this one must go to scales smaller than 
the free-streaming scale of dark matter particles. For 
example, the comoving free-streaming scale of a typical 
neutralino dark matter particle is of the order of parsec. 
We do not study this system on larger cosmologically 
relevant scales. It is, however, unlikely that the system 
will still able to drive the accelerated expansion of the 
universe since the generic AZK instability is intimately 
related to the adiabatic sound speed of the fluid fl^ ]. 

The result above assumes that the scalar field is much 
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heavier than the expansion rate of the universe. This 
constraint is easy to satisfy and the large mass makes the 
model more robust to radiative corrections (for example, 
see (HI). Secondly, the calculation is only valid for modes 
which have a wavelength much larger than l/m x ; for 
shorter wavelengths, we cannot assume that the scalar 
field relaxes to the minimum quickly enough. 



AVOIDING THE AZK- INSTABILITY 

To avoid this instability, we look at more general cou- 
plings. 

Consider now a model where the total energy is 

E = J d 3 x [V (x) + mn M + Xg{ X )n^] , (7) 

and we choose A > without loss of generality. 

Again we assume that the scalar field tracks the mini- 
mum of the potential and hence we have, 

K(xo) + A.g'( X oKi = 0,(8) 

K( X o) + A 5 "(xo«) (§^) + Xg\ XQ )nn^ = .(9) 

Following our earlier calculation, we find 



af- 1 [^^mVi ["V(Xo)»S] 2 



+Xn(n - l)g(xo)n 



m I > 



(10) 



Therefore, the instability is avoided if 



n 2 AV [5 ' (X ° )r 



l M" 2 

mt 



+ n(n-l)A 5 ( Xo )>0. (11) 



We note that the first term is always negative and 
gets large with hm unless g'(xo) decreases fast enough. 
Looking at the second term we note that any value of 
< n < 1 is unstable independent of the form of g(x) ex- 
cept for the requirement that g{xa) > which is required 
anyway for the potential to be bounded from below. 

A robust way to avoid the instability is to choose n < 0, 
which makes the second term positive. This is, of course, 
not sufficient to guarantee the inequality in Eq. QT] We 
need the magnitude of the second term to be larger than 
that of the first. This is easy to arrange. We again look 
at changes to the potential as we vary \ about xo- If the 
potential is not fine-tuned to give rise to cancellations be- 
tween terms in the Taylor expansion, then nXg'n^dx < V 
and also m 2 ,(6x) 2 /2 < V. Putting these two expressions 
together yields 2n 2 X 2 (g') 2 n 2 ^ /m 2 , < V ~ Xgn^. Hence 
we see that it is natural, if n < 0, for the inequality in 
Eq. [11] to be satisfied. 



It is also possible to avoid the instability by choosing 
n > 1. However, this region of model space will be heav- 
ily constrained by observations. In situations where the 
matter density gets large, i.e., in collapsed structures, the 
last term in the potential dominates. It would make the 
dark energy density in galaxies large, change structure 
formation and clustering properties of dark matter ha- 
los. Therefore, these kinds of models would be tightly 
constrained. In order for these models to be viable, A 
would have to be small and the model would essentially 
be the same as that with two decoupled fluids. 

Thus the requirement of AZK-stability and observa- 
tional constraints naturally lead us to consider models 
where n < 0. We now look at observational consequences 
of such a coupling. 



AZK-STABILITY AND SUPER-ACCELERATION 

The coupling term above with n < introduces a 
very interesting effect: this model has super-acceleration. 
That is, observations will seem to show a phase with dark 
energy equation of state less than -1. 

To see this, we first note that the observational quan- 
tity that is important is the pressure. We will fit to the 
observations a model with matter that scales with the 
expansion as 1/a 3 , and dark energy with some equation 
of state wde- Note that adding or removing a component 
of energy density that scales as 1/a 3 does not change the 
pressure of the fluid. Hence very generally Ptot = Pde- 
P to t is defined by the equation V — —3H(V + Ptot) from 
which we find P tot = -V ( X o) + Xg( Xo )n^(n-i). We set 
the equation of state wde = Ptot/iV — mnyi) and find, 



1"DE 



-1 + 



^Ag(xo)»M 
^o(xo) + Ag(xKi 



(12) 



Now since n < 0, the second term is actually negative, 
and we have wde < —1 i-e. super-acceleration. 

We emphasize that this super-acceleration is not ac- 
companied by any of the problems normally associated 
with theories with equation of state less than -1. There 
is no acausal behavior, and there are no ghosts. This is 
because the super-acceleration in our model results from 
an interaction which is ignored in the fitting of theory to 
observations. If we fit our observations using a canonical 
matter density term and dark energy, then the interac- 
tion has the effect of making the the effective equation of 
state for dark energy more negative. 



SOUND SPEED 

Here we present an alternative derivation of the insta- 
bility in terms of the sound speed of the combined fluid. 
A negative sound speed squared would signal instability. 
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On length scales much larger than m x 1 , the evolution 
of the system is adiabatic and hence the sound speed is 



"tot 



(13) 



The adiabatic sound speed in this theory can then be 
expressed as 



"-mi9wde/(9"m + wde(1 + wde) 
1 + wde + mn M / (V - mn M ) 



M 



d 2 V(xo,n M ) 



( dxo 



x \dn M 
n M dwt t/dn M + u>tot(l + w tot) 



1 + Wtot 



(14) 
(15) 
(16) 



where Wtot = Ptot/V is the equation of state of the total 
fluid. 

For a universe with an accelerating expansion u> to t < 
— 1/3. For a wide class of models with iy to t < and either 
the nMw'tot term sub-dominant or negative, we have c 2 < 
and the system is unstable. This is just the AZK- 
instability. 

Lets now look in more detail at Eq. [T3J First, consider 
the case where u>de > — 1: the denominator is positive 



and if the 
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term is sub-dominant or negative, then 



AZK-instability sets in. It is clear that this instability 
may not be present in models with wde < — 1. We also 
note that this instability will likely set in well before the 
current epoch because at early times Um/ {V — Um) 3> 1- 
For this case where wde(1 + wde) > 0, the sign and mag- 
nitude of the n-M^DE term is important. In particular, 
the requirement that the nMWn R term is sub-dominant 
may not be trivial to obtain 61 1. 

While the above derivation shows us how the instabil- 
ity arises, it does not provide us with an intuitive un- 
derstanding of what happens to the matter. In order to 
better understand that we look at the Boltzmann equa- 
tion for the matter coupled to a scalar field. The scalar 
field gives the matter a mass term that can vary spatially 
and temporally. Following AZK we write down the 
Boltzmann equation for matter neglecting gravity and 
hence only valid on small scales. These are the scales of 
interest since we have assumed m x 3> H. We write down 
the first order perturbations to this equation and expand 
the perturbations in plane wave modes. Denoting the 
effective mass of the matter particle by M(x) we find, 

«a/(p,k) - ( 7 M)- 1 P -w(p,k) (17) 

+ 7 - 1 < 5M(k)k- V p /(p) -0. (18) 

We then find the perturbation to the matter density 
<5n,M(k) using the above equation. In the limit that mat- 
ter is non-relativistic, the resulting equation has a simple 
form. We find that the variation in effective mass of the 
particle is given by SM(k) = (M/riM)c 2 dn,M(k) where we 



have defined c s = cu/k, the sound speed of matter. The 
above equation is valid for perturbations SM on all scales 
at which our assumptions hold. As pointed out in [lij ]. 
there is no scale in the equation for c 2 because we are 
studying scales where it is correct to assume that the 
scalar field adjusts to changes in the matter density, and 
gravity is unimportant. 

We now turn to the fluid description and write M = 
V(Xo, ?im)/<9«m- Using Eq.[4]for dxo/driM, one may then 
obtain perturbations in M as SM — [M / ' nyijc^Snyi where 
c 2 is given by Eq. 1151 In the framework of a scalar degree 
of freedom coupled to matter, both descriptions must be 
valid and hence we find that c 2 = c 2 . The instability may 
therefore be analyzed in terms of c 2 . All of our analyses 
in earlier sections go through if we work with c 2 and we 
conclude that models with super-acceleration provide a 
generic way to avoid the AZK instability. 



CONCLUSIONS 

In this paper, we have explored the possibility that 
dark energy may interact with matter. Such a hypothesis 
is natural if the explanation for dark energy requires ex- 
tra scalar degrees of freedom. Unfortunately, as we have 
shown here, these models suffer from a generic instabil- 
ity when the mass of the scalar field is very large. We 
have verified that this instability is also present in scalar- 
tensor theories where the scalar plays the role of dark 
energy, and also in models with multiple scalar fields. 

We then looked for models where this instability could 
be avoided, and found a large class of such models. Most 
interestingly, we found that in these models, the apparent 
equation of state of the dark energy density is generically 
smaller than -1. This super- acceleration is a result of the 
fact that we fit observations with models that have non- 
interacting matter and dark energy fluids. 

There is a theoretical prejudice against models of 
wde < — 1 due to their apparent theoretical problems. 
The observational data certainly do not disfavor u>de < 
— 1. Indeed a large region of the parameter space al- 
lowed by SNIa observations corresponds to a constant 
wde < — 1. Here we have shown that stable models with 
wde < — 1 may be constructed without encountering 
ghosts or acausal behavior. These models are no more 
fine-tuned than quintessence models. Thus theoretical 
bias against wde < — 1 should be treated with circum- 
spection, and not be given any weight when interpreting 
observational data. 
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